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We consider an asymptotically free vectorial SU(Afc) gauge theory with Nf fermions in the funda- 
mental representation and analyze higher-loop contributions to the evolution of the theory from the 
ultraviolet to the infrared in the limit where Nc ^ oo and Nf — > oo with r = Nf /Nc a fixed, finite 
constant. We focus on the case where the n-loop beta function has an infrared zero, at ^ = £,ir^„i, 
where ^ — aNc. We give results on ^iR^„e, the anomalous dimension of the fermion bilinear eval- 
uated at (,iR^ne, denoted ^iR^nt, and certain structural properties of the beta function, /3g. The 
approach to this limit is investigated, and it is shown that the leading correction terms are strongly 
suppressed, by the factor 1/Nc ■ This provides an understanding of a type of approximate universal- 
ity in calculations for moderate values of Nc and Nf, namely that aiR^„tNc, 'yiR,ni, and structural 
properties of the beta function are similar in theories with different values of A^c and Nf provided 
that they have similar values of Nf /Nc- We give results up to four loops for nonsupersymmetric 
theories and up to three loops for supersymmetric theories. 

PACS numbers: 11.15.-q,11.10.Hi,11.15.Pg 
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I. INTRODUCTION 

The evolution of an asymptotically free gauge theory 
from the ultraviolet (UV) to the infrared (IR) is of funda- 
mental field-theoretic interest. The UV to IR evolution 
of the gauge coupling g{^) as a function of the Euclidean 
momentum scale, /x, is determined by the j3 function [ij 



B = — 



(1.1) 



where t = In/i and a{p) = gip)"^ /{An). Here we con- 
sider this evolution for a vectorial gauge theory with 
gauge group G — SU(A^c) and Nf massless fermions ^pj, 
j — 1, Nf, transforming according to the fundamental 
representation of G [H . We also point out some contrasts 
with results for fermions in higher-dimensional represen- 
tations. We focus on the case where the n-loop /? function 
has an infrared zero, at a value a — am^ne- The condi- 
tion of asymptotic freedom requires that Nf be bounded 
above by a value where the one-loop coefficient in 

/? vanishes W\- For large enough Nf (less than Nf^tiz), 
the two-loop /3 function has an infrared zero at a certain 
value of a, denoted aiR,2e [1, The desire to under- 
stand better both the behavior of the running coupling in 
quantum chromodynamics (QCD) and the properties of 
an IR zero that occurs for sufficiently large Nf have mo- 
tivated calculations of higher-loop terms in /] [1, 0] and 
higher-loop corrections to the two-loop result for the IR 
0. In 1,113, ons of the IR zero in /?, 



zero 

and the associated anomalous dimension of the (gauge- 
invariant) fermion bilinear, 7™, were done to four- loop 
order for an asymptotically free vectorial gauge theory 
with gauge group G and Nf fermions in an arbitrary 
representation R, with explicit results for R equal to the 
fundamental, adjoint, and symmetric and antisymmet- 
ric rank-2 tensor representations. Further generalizations 
and results on higher-loop structural properties of the f3 
function were given in [13 |. 



An interesting and important property that one no- 
tices in these calculations for an SU(A'^c) theory with Nf 
massless fermions in the fundamental representation is 
that the values of the n-loop 7™, evaluated at aiR^ne, 
denoted jiR,ni, and of the product aiB.,niNc, are simi- 
lar for theories with different values of Nc and Nf, pro- 
vided that these theories have similar values of the ratio 
Nf/Nc- Indeed, the computations in |12i] show that this 
is also true for other structural quantities describing the 
UV to IR evolution, including the derivative of the n- 
loop beta function, d/3ne/da evaluated at am^ne and the 
products ar,i,neNc and {Pne)minNc, where am,n£ denotes 
the value of a where f3ni is a minimum, and (/3n^)mm is 
the value of (3ne at this minimum. These observations 
show that there is an underlying approximate universal- 
ity in the form of the quantities that control the UV to 
IR evolution of these theories. This motivates a more 
detailed study to elucidate this phenomenon. We carry 
out this study in the present work. 

For this purpose, we analyze these theories in the 't 
Hooft - Veneziano limit jisl . [3| 

Nc^ 00 , Nf 00 



with r ; 



Nc 



K and ^(/i) = a{p)Nc = A(/i) 



(1.2) 



where «; is a constant and A is a function depending only 
on fi. We will use the symbol lim^jvw for this hmit, 
where "LNN" stands for "large Nc and Nf" (with the 
constraints in Eq. (II. 2p imposed). The reasons for these 
two constraints in Eq. (|1.2p (that the ratio Nf/Nc and 
the product a{fj,)Nc are fixed and finite) are that these 
constitute the necessary and sufficient conditions in order 
that (i) fermions give nonvanishing contributions to the /3 
function, anomalous dimension 7^, and other quantities, 
and (ii) scattering amplitudes remain finite, in the limit 



2 



TVc — > oo, respectively. More generally, if the fermions 
are nonsinglets under other gauge groups with squared 
couplings Ui, one also_jequires that the products UiNc 



be finite as Nc ^ oo [15] 



As we will show in detail below, a study of the LNN 
limit p. 211 and the approach to it provides an explanation 
of the approximate universality in the UV to IR evolu- 
tion of theories with different values of Nc and Nf but 
the same, or similar, values of r that is exhibited in ex- 
plicit calculations (with appropriate scalings understood 
for certain quantities, such as multiplying aiR^ni by Nc). 
A crucial property of the LNN limit is that it reduces the 
number of variables on which the UV to IR evolution de- 
pends. Thus, for finite and Nf, this evolution and the 
j3 function that describes it, depend on three variables: 
a (and thus, parametrically, /Lt), and TV/, while in the 
LNN limit, they only depend on the two variables a(/i) 
and r. 

Since the rational numbers Q are dense in the real 
numbers R, it follows that in the LNN limit, one can 
choose values of N^. and Nf so that the rational number 
r is arbitrarily close to any non-negative real number. 
Therefore, henceforth, to arbitrarily good accuracy, we 
may simply treat r as a real number, and we will do so. 
As is well-known, in the Nc — oo limit and also in the 
LNN limit, the gauge group SU(A^c) is effectively equiv- 
alent to \]{Nc). The use of a large-iV limit, where N is 
the number of components in a spin or field, has been 
valuable in the past partly because it allowed one to ob- 
tain exact results for statistical mechanical models (T6j 
and quantum field theories [l^- 15 1, [13, [3. Our pur- 



pose in using it here is somewhat different, namely to 
gain further insight into the above-mentioned approxi- 
mate universality that is exhibited by calculations of the 
UV to IR evolution of theories with different Nc and Nf 
but equal or similar values of r. 

As part of our analysis, we will briefly contrast the 
properties of theories with fermions in the fundamental 
representation with properties of theories with fermions 
in higher-dimensional representations. In the case where 
fermions are in a two-index representation (including the 
adjoint, and symmetric and antisymmetric rank-2 tensor 
representations), the condition that is necessary and suf- 
ficient to construct a finite Nc ^ oo limit is to set Nf 
equal to a (non-negative, integer) constant. This is a re- 
sult of the fact that the quadratic Casimir invariants Tf 
and Cf [l^ that enter into the coefficients of the beta 
function grow like (a constant times) Nc as Nc —?' oo. 
If the fermions are in a representation involving three 
or more indices, then generically for fixed finite Nf, the 
fermion contributions dominate over the gauge field con- 
tributions to the P function by powers of Nc as Nc ^ oo. 
For example, for the symmetric rank- 3 tensor represen- 
tation, Tf ^ N^ for large Nc, so that the fermion con- 
tribution to the leading (3 function coefficient dominates 
over the gauge-field contribution, which is ~ Nc, spoiling 
the asymptotic freedom of the theory. Hence, aside from 
our primary focus on the case of fermions in the funda- 



mental representation, we will restrict our discussion of 
other representations to the adjoint, and symmetric and 
antisymmetric rank-2 tensors. 



By taking r near to its maximum value allowed by 
asymptotic freedom, one can arrange that the zero of /3 
occurs at an arbitrarily small value of and one may 
conclude that in the infrared the theory is in a decon- 
fined non-Abelian Coulomb phase without any sponta- 
neous chiral symmetry breaking. In this case, the IR 
zero of /3 at is an exact fixed point of the renormal- 
ization group for the theory. In contrast, as r decreases, 
^iji increases, and studies with finite Nc and Nf lead 
to the conclusion that for Nf less than a critical value, 
Nf^cr, as fj. decreases though a value denoted A, the in- 
teraction strength a(/i) exceeds a critical value, acr, to 
produce spontaneous chiral symmetry breaking and asso- 
ciated dynamical mass generation for the fermions. For 
a given Nc, the theory may thus be considered to un- 
dergo a (zero-temperature) chiral phase transition as Nf 
passes through this value, Nf^cr [2Q[, and there has been 
an intensive research program using lattice gauge simu- 
lations to determine A^/.cr for values such as Nc — 3 and 
Nc = 2 [2lj. Correspondingly, in the LNN limit consid- 
ered here, the theory undergoes a chiral phase transition 
as r passes through rcr, where rcr = Nf^cr/Nc, with the 
UV to IR evolution leading to a chirally symmetric phase 
for r > rcr and a phase with spontaneous chiral symme- 
try breaking for r < rcr- If r < rcr, then in the effective 
low-energy field theory below A, one integrates out the 
fermions (which have dynamically generated masses of 
order A), and the /3 function changes to become that of 
a pure non-Abelian gauge theory, which does not have a 
(perturbative) zero. In this case, S,ir is only an approxi- 
mate fixed point. 



This paper is organized as follows. In Section |TI] we 
define an appropriately scaled beta function, called /3^, 
that is a finite function of ^ in the LNN limit and in 
this section, and in Sections IIIII and HVl we investigate its 
structure up to four-loop order. Our results include an 
analysis of the behavior of the coefficients as functions 
of r, the LNN limits for the n-loop IR zero, £,iR,ne, the 
value of ^ where is a minimum, the value of at 
this minimum, and the derivative d/J^/d^ evaluated at 
£.iR,n£- III Section|V]we carry out a similar analysis of the 
coefficients in the anomalous dimension of the fermion 
bilinear, 7^, and its value at ^/_r,„£, again up to four- 
loop order. In Section IVII we calculate correction terms 
to the LNN limits for various quantities and give a general 
analytic explanation for the rapidity with which this limit 
is approached, namely that these correction terms are 
strongly suppressed, by the factor 1/N^. Section IVlIII 
is devoted to a corresponding study of the LNN limit 
of a supersymmetric gauge theory. Our conclusions are 
contained in a final section. 
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II. I3 FUNCTION AND SOME GENERAL 
PROPERTIES IN THE LNN LIMIT 

A. General 

In this section we analyze the /3 function in the 
Umit (|1.2|) . It win be convenient to define a(/i) = 
5(m)V(16^') = a(M)/(47r) and 



47r 



(2.1) 



(The argument /i will often be suppressed in the nota- 
tion.) In terms of a, or equivalently, a, the beta function 
has the series expansion 



/3 = /3a = 



oo 



£=1 



i=i 



(2.2) 



Thus, similarly to the relation between bg and bf 



bi 



(47r)^ 



(2.6) 



As with Eq. (|2.5p . it is understood here and below that 
all expressions have been evaluated in the LNN limit 
p.2p . The function, calculated to n-loop {ni) order, 
is denoted by fi^^ni and is given by Eq. \2A\ with oo 
replaced by n as the upper limit on the sum over £. 

To analyze the zeros of /S^^nf, aside from the double 
zero at f = X = 0, we extract an overall factor of —2^^ 
and calculate the zeros of the reduced (r) polynomial 



2e 



1 

47r 



1 X 



(2.7) 



where i denotes the loop order and b^ = 6£/(47r)^. Thus, 
the n-loop beta function is given by Eq. ()2.2p with oo 
replaced by n as the upper limit on the summation over 
loop order, The coefficients bi for ^ = 1, 2 are in- 
dependent of the scheme used for the regularization and 
renormalization of the theory and were calculated in [3,] 
and [il. The bi with £ > 3 are scheme-dependent and 
have been calculated up to {i = 4)^loop order d, Bin 
the modified minimal subtraction [22] (MS) scheme [2a |. 
The usefulness of the AIS scheme has been demonstrated, 
e.g., by the fact that inclusion of three-loop and four-loop 
corrections in the running of as(/i) in QCD significantly 
improves the fit to experimental data [2J|. The scheme- 
dependence of the higher-loop IR zero of the beta func- 
tion was recently studied in |25| . 

For our present analysis of the theory in the LNN limit, 
the first step is to construct a beta function that has a 
finite, nontrivial LNN limit. We do this by multiplying 
both sides of ([^ by and then taking the LNN limit. 
The result is a function of ^ and can be expressed as 



13.= ^^ lim 13, 



Nr. 



(2.3) 



This function has the expansion 

oo oo 



As is clear from Eq. (|2.7|) . the zeros of /3^^ne away from 
the origin depend only on rt — 1 ratios of coefficients, 
which can be taken as bi/bn tor £ — 1, n — 1. Although 
Eq. (|2.7I) is an algebraic equation of degree n — 1, with 
n — 1 roots, only one of these is physically relevant as the 
IR zero of P^^ne- We denote this as ^/_r,„£ — 47ra;^^ 



B. Behavior of Coefficients in /3 as Functions of r 

From the expressions for bi and 62 [1, 0| i have 



and 



61 = -(11 -2r) 



62 = -(34-13r) . 



(2.8) 



(2.9) 



where 



In the MS scheme, from the expression for 63 [6|, we 
obtain 



63 = — (2857- 1709r-|- 112r2) 
54 

= 52.9074- 31. 6481r + 2.07407r2 (2.10) 



lim — 7- 

LNN 



lim — 7- 

LNN N^. 



(2.5) 



and from 64 |7|, we obtain 



f 150473 /485513\ /8654n , /130\ ,4, . 2nw.n 

64 = ( ]r+{ r^+l r^-l-- Il-5r-f21r2 C 3) 

^ 486 V 1944/ V 243 / V243/ 9^ '^^ ' 

= 315.492 - 252.421 r -I- 46.832 r^-f- 0.534979 , (2.11) 
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to the indicated numerical floating-point accuracy, where 
C(s) = is the Riemann ( function. For some 

purposes it is more convenient to deal with the bi, since 
they are free of factors of 47r, while for numerical purposes 
it is often more convenient to use the since the range 
of values of bi as functions of £ is somewhat smaller than 
the range for the bi. In Table U we list values of bi for 
2 < £ < 4 as functions of r in the interval < r < rf,iz. 

In 0, Nf, btz was defined as the value or set of values 
of Nf where b^ — and thus, for our present analysis, we 
define 



= hm — 

LNN Nr 



From Eq. (|2.8|) . we have 



11 

y 



(2.12) 



(2.13) 



As is evident from Eqs. (12. 8p and (|2.9I) . the coefficients 
bi and 62 are both monotonically (and linearly) decreas- 
ing functions of r. The coefficient 61 decreases from 11/3 
to as r increases from to ri,iz — 11/2. We require 
that the theory be asymptotically free, i.e.. 



r < niz 



11 



(2.14) 



The coefficient 62 decreases from 34/3 at r = and 
passes through zero to negative values as r increases 
through the value 



'rb2z = 777 



34 
13 



As r Thiz, &2 reaches the value 
^ 25 

62 = — — at r = niz ■ 



(2.15) 



(2.16) 



Therefore, in the LNN limit, the interval in r, denoted Ir, 
where the two-loop /3j function has an IR zero, is given 
by 

Ir ■■ n2z <r< niz , < ^ < y (2-17) 

(i.e., 2.615 < r < 5.500). With r € Ir, we will, cor- 
respondingly, focus on the UV to IR evolution in the 
interval 

h- < ^(Ai) < ^iRM . (2.18) 
The coefficient 63 vanishes at two values of r, denoted 
1709 - 57\/5()5 



''632 1 



224 



= 1.911 



and 



1709 + 57v/505 
rb3z2 = ^ = 13.348 , 



(2.19) 



(2.20) 



where here and below, the floating-point values are given 
to the indicated accuracy. This coeflticient 63 is mono- 
tonically decreasing in the interval < r < Vbiz, decreas- 
ing from 63 = 2857/54 = 52.907 at r = and passing 
through zero to negative values as r increases through 
fbazi in Eq. (|2.19p . As r increases from nszi, &3 contin- 
ues to decrease and passes through the value 



5299 
338 



= -15.6775 at r : 



rb2z 



34 
13 



(2.21) 



at the lower end of the interval Ir- As r increases 
throughout the interval Ir, b^ decreases further, and as r 
increases to its maximum, rbiz, at the upper end of this 
interval, 62 reaches the value 



63 
Since 

and 



701 



-58.417 at r = 7- 



rb3zl < n2z 



11 



biz 



rb3z2 > niz , 

it follows that in the MS scheme, 
63 < y reir 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



Given this result and the fact that the quantity —54^3 
will appear in later formulas, it will be convenient to 
denote 

Dse = -5463 = -2857 + 1709r - 112r^ , (2.26) 

which is positive for r Ir- 

For completeness, we note that 63 reaches a minimum 
at r = 1709/224 = 7.629, and, for larger r, it increases, 
passing through zero again at the value rf,3z2 in Eq. 
(|2.20p . Since these values of r lie above rbiz, they are 
not of direct interest for our present study. 

As r increases through the range < r < niz, the 
coefficient 64 in the M S scheme decreases from the value 

^4^^ + ^- 315.492 at . = 0, (2.27) 



486 



9 



(i.e., 64 = 1.265 X 10 ^), passes through zero with nega- 
tive slope at 



rw. 1 = 2.040 



(2.28) 



reaches a minimum of —14.831 at r = 2.581, and then 
increases. At the lower end of the interval Ir, at n2zj 

f 550009 31900 C(3) 

64 = — — -I -f!-^ = -14.770 at r = r^z ■ 



6084 



507 



(2.29) 
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As r increases in the interval 7^, 64 passes through zero 
again, at 

= 3.119 , (2.30) 
this time with positive slope, and attains the value 
14731 

64 = -j^ + 275C(3) = 432.864 at r = (2.31) 

(i.e., 64 = 1.736 X 10~^) at the upper end of the interval 
Ir- Some special values, expressed in terms of the hg 
coefficients, are listed in Table HIl 

Concerning the sign of 64 in the range r > 0, 

64 > for < r < TbAzA , 

64 < for rb4z,i <r < rbiza , 

64 > for rMz,i <r < niz . (2.32) 

That is, numerically, S4 > if < r < 2.040 or r > 
3.119, and 64 < if 2.040 < r < 3.119. The zero of 64 
at rf,42,i = 2.040 lies below the lower end of the interval 
Ir (at ri,2z = 2.615), while the zero at rb4z,2 = 3.119 lies 
in the interior of the interval Ir- Hence, restricting to 

For r e Ir , 64 < if rb2z <r < tmz,! 



64 > if niz,! <r < niz , 

(2.33) 

i.e., numerically, for r e 7^, 64 < if 2.615 < r < 3.119 
and 64 > if 3.119 < r < 5.500. Since 64 is a cubic poly- 
nomial in r, there is a third value of r where it vanishes, 
but this is at the negative, and hence unphysical, value 
r = —92.699 and hence is of no direct relevance here. 
Although our analysis here presumes the LNN limit, a 
remark is in order for finite Nc and Nf. The interval 
where 64 is negative for r e is not present for suffi- 
ciently small Nc and Nf. This is evident from the ex- 
plicit 64 values listed in Table I of our Ref. [9] for Nc = 2 
and Nc — 4. This interval of negative 64 values is present 
for Nc > 4. 



III. IR ZERO OF /3 

Combining the results from the previous section, we 
exhibit the explicit four-loop function. For this pur- 
pose, it is simplest to use the x variable defined in Eq. 
([2T|) . We have 



= -Sttx^ 



11 - 2r 



34 - 13r 



2857- 1709r+ 112r2 
54 
4 



r 150473 /485513\ /8654\ , /130\ .4, ^l ^ 4^ 



r 



(3.1) 



A. Two-Loop Level 



At the two- loop level, if r e , then has an IR zero 



at 



too large for this perturbative calculation to be applica- 
ble. It will be useful here and below to give values of 
various quantities at an illustrative value of r. 



21 



47r(ll - 2r) 
13r - 34 



(3.2) 



Since this is obtained from a perturbative calculation, it 
is only reliable if S,ir^21 is not too large. As r — >■ r^iz, 
^iR,2e — 0, and hence in the upper end of the interval 
Ir, one may plausibly expect that this two-loop expres- 
sion becomes a progressively more and more accurate ap- 
proximation to the IR zero of the exact /3j function. As 
r \ 34/13 at the lower end of the interval Ir, £,iR,2i grows 



B. Three-Loop Level 

At the three-loop level, the IR zero of /3{ is given by the 
physical (smallest positive) root of the quadratic equa- 
tion 

I3i.2e,r = h+ 62^ + ^3^' = . (3.3) 
This equation has, formally, two solutions, namely 



263^ 



-( - 62 ± 46163 



(3.4) 
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Since we have shown that 63 < for r G Ir for a general 
scheme that preserves the existence of the IR zero in the 
(scheme-independent) /3j,2^ at the three-loop level, we 
can rewrite (|3.4p as 



-^(-|62|TvW+4^ 



(3.5) 



As is evident from Eq. (|3.5p , only the root corresponding 
to the lower sign choice in Eq. p.Sp is positive and hence 
physical. We denote it as 



1 



216, 



(3.6) 



By the same type of proof as was given in [^[T3l, for the 
relevant interval r € Ir where the scheme-independent 
two-loop l3^^2£ function has an IR zero, we find that at 
the three-loop level 



IR,3£ 



IR,2£ 



for r £ Ir 



(3.7) 



with equality only at r = r&iz, where ^7^,3^ = ^ma^ = 0. 
In we pointed out that the corresponding inequality 
OiiR,u < 0iiR,2i applies more generally than just in the 
MS scheme, and the same is true of the inequality (|3.7p . 
We recall the reasoning for this. Since the existence of 
an IR zero in the two- loop j3 function, P^^2e, is a scheme- 
independent property of the theory, a reasonable scheme 
should maintain the existence of this IR zero (albeit with 
a shifted value) at higher-loop order. Now in order for a 
scheme to maintain this zero, a necessary and sufficient 
condition is that 63 ^ 46163 > 0, so that the square root 
in Eq. (|3.4p is real. But the lower end of the interval 
Ir is defined by the condition that 62 — )■ as r \ ri)2z- 
Given that r S /r so /3^^2e has an IR zero, this means 
that a reasonable scheme, which preserves the existence 
of this zero at the three-loop level, should have 63 < 
for r G Ir- From this, by the same type of proof as was 
given in for this class of schemes, the inequality (13. 7p 
follows. 

Substituting the relevant expressions for the 6^ in 
we have, in the MS scheme, the explicit result 



^IR,M 



127r[-3(13r 



(3.8) 



where D^g was defined above in Eq. (j2.26p . and it is 
convenient to define the shorthand notation 



52450 + 41070r - 7779r^ + 448r'^ 



(3.9) 



The polynomial C^g has only one real zero, at r = 
1.86532 (to the indicated accuracy) and is positive for 
r > 1.86532, and hence for all r S Ir- The polynomial 
— 3(13r — 34) vanishes at the lower end of the interval Ir 
and is negative for r £ Ir, but it is smaller than \fCu, so 
^iRM > for r e /r, as is necessary for it to be physical. 

In Table mil we list numerical values of S,ir^21 and ^iB.,3e 
for r Cz Ir- As is evident in this table, (,ir^21 and £,irm 



decrease monotonically as a function of r throughout this 
interval Ir (as does to be discussed below). At the 

lower end of this interval. 



C/fl,3£ = 207r 



26 
5299 



= 4.401 



and at the upper end, 
ilR 



34 

at r = — 
13 



11 



-> a.s r /• niz Y 



(3.10) 



(3.11) 



The ratio ^iR,zi/S,iR,2i increases monotonically from as 
r increases from the value r = 34/13 at the lower end of 
the interval Ir , and this ratio approaches 1 from below as 
r approaches the upper end of the interval at r = 11/2. 
It is useful here and below to give illustrative values of 
various quantities and ratios at an illustrative value of r. 
For this purpose, we choose an r approximately in the 
middle of the Ir, namely r = 4. We have 



iiR,: 



2tt 



2.0944 



and 



4^(-2 + V22) 

?/i?,3£|r=4 — 27 ~ 



2522 



so that 



^iR,3e 



£,IR,2l 



2(-2 + V22 ) 



= 0.5979 



(3.12) 



(3.13) 



(3.14) 



This ratio provides an illustrative measure of the decrease 
in the value of the IR zero of /3 when one calculates it at 
three-loop order, as compared with two-loop order. 



C. Four-Loop Level 

At the four-loop level, the IR zero of is the (smallest 
positive) root of the cubic equation 



61 + 62^ + 63^ + he = 



(3.15) 



Now 62 < for r E Ir, and we recall our discussion 
above, that 63 < for r E Ir in the A/5 scheme and 
other schemes that maintain the existence of the IR zero 
in (3^.21 at the three-loop level. We can therefore write 
Eq. (1X151) as 



6i-|62|e-|63|C' + ^4C' = 



(3.16) 



For r E Ir, Eq. p.l5p . or equivalently, (I3.16p . has three 
real roots, and from these we determine the relevant 
(smallest, positive) one as £,ir,m- We list values oi£,iRAi 
in Table Uni 
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D. Shift of IR Zero From n-Loop to (n + l)-Loop 
Level 



In [121, ^ general result was derived concerning the 
sign of the shift of the IR zero of /3 going from the n- 
loop level to the {n + l)-loop level. Provided the scheme 
has the property that bi with € > 3 are such as to 
maintain the existence of the zero in the two-loop /3 
function, then a/^ („+i)^ > am^ni if if bn+i > and 
ctiR,{n+i)e < oiiB.,ni if if ^n+i < 0. Thc same proof can 
be applied here to deduce that, provided that the scheme 
has the property that bt with ^ > 3 are such as to main- 
tain the existence of the zero in the two-loop function, 
then 

ilR,{n+l)l > S,IR,ni if hn+1 > 0, 
^IR,{n+l)e < aiR,ni if bn+1 < . 

(3.17) 

We may apply this inequality for the comparisons of 
S.iR,3e with C/-R,2£ and ^ir^h with ^ir^si- Since 63 < 
for r ^ Ir, this result provides another way of deduc- 
ing the inequality (13.71) for the two-loop versus three- 
loop comparison. Applying the general inequality for the 
three-loop versus four-loop comparison, we infer that 

^IRM < ^iR,3e if 2.615 < r < 3.119, so 64 < 0, 

^iRAi > ^IRM if 3.119 <r< 5.500, so 64 > . 

(3.18) 

These inequalities are evident in Table Hill For example, 
at r = 3.0, £,iR,u/£,iR,u = 0.976, while for r — 5.0, 
^iRAtl ^IR,M — 1-02. One sees that the magnitude of the 
fractional difference 



le 



ilR,M\ 



S.IR,U 



(3.19) 



is reasonably small. This is in agreement with one's gen- 
eral expectation that if a perturbative calculation is reli- 
able, then as one calculates this quantity to progressively 
higher-loop order, the magnitudes of the fractional dif- 
ferences between the values at the n'th and {n + l)'th 
orders should decrease. 



E. Summary of Results on IR Zero of 

We summarize our findings concerning ^ib.m fol" 
lows. As one goes from the (scheme-independent) two- 
loop level to the three-loop level, the value of the IR zero 
of /3j decreases. For r in the lower part of the interval 
Ir where the two-loop function has an IR zero, this 
reduction in the value of the IR zero is rather substan- 
tial. For example, for r = 3.0, near the lower end of the 



interval Ir, £,ir.m/£,ir,21 = 0.234, while for r = 5.0, near 
the upper end of I^, £,1 r,u / £,1 R,2i — 0.873. Going from 
three-loop to four-loop order, the change in the value of 
the IR zero is smaller in magnitude and can be of either 
sign, depending on the value of r G . In general, both 
iiR,u and iiR^t are smaller than C/_r,2£- 



IV. SOME STRUCTURAL PROPERTIES OF /^^ 

For theories which exhibit an IR zero, S.iR^2e, in the 
two- loop beta function, /3^,2£, there are several structural 
properties of interest in addition to higher-loop values of 
this IR zero. These include 

• the value of 1^ at which l3^^ni reaches a minimum in 
the interval /j, denoted £m,ni, where the subscript 
m denotes minimum 

• the value of „£ at this minimum, denoted 

(/?«,«£) mm 

• the derivative of P^^nt at £,ir^„i, denoted 



iJR,ne 



d£ 



(4.1) 



Note that because £ — aNc and (3^ — limijvAr PaNc, the 
factor of Nc divides out in the derivative df3^/d£, so that 



-T^ = hm — — 

dS, LNN da 



and 



In particular 



d£ 



— lim 



dj^a.r, 



LNN da 



d£ 



lim 



i=^IR.,^t LNN da 



(4.2) 



(4.3) 



(4.4) 



As was discussed in [12], higher-loop calculations of the 



derivative 



are of interest because this en- 



da |a=a/H., 

ters into estimates of a dilaton mass in a quasiconfor- 
mal gauge theory. In turn, this also provides one mo- 
tivation for studying the LNN limit of this derivative. 



A. Position of Minimum in /?{,n< 

Concerning the position of the minimum in for r S 
Ir, we calculate that at the two- loop level. 



£m, 



2i 



87r(ll - 2r) 
3(13r- 34) 



This satisfies 



im,2i — -^£iR,: 



(4.5) 



(4.6) 
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At the three-loop level in the M S scheme, we find 
37r[-9(13r-34) + 



Dm 



(4.7) 



where we define the shorthand notation 

Eu = -409196 + 320604r - eOTllr^ + 3584r^ 
From Eqs. ((3^ and (|4Jl) . we find 

U,U -9(13r-34) + V^ 



i 



4[-3(13r-34) + VC^] 



(4.8) 



(4.9) 



This may be compared with the corresponding ratio of 
two-loop quantities ^m,2i/^iR.2i — 2/3. In contrast to the 
latter ratio, which is a constant, independent of r S Ir, 
the ratio (14.91) is a monotonically decreasing function of 
r £ Ir- At the lower end of this interval, 



Um 1 , 34 
= -= at r = r&2. = -r 



(4.10) 



As r approaches the upper end of the Ir at r^iz = 11/2, 
the ratio ()4.9p approaches the limit 



lim 



(4.11) 



Note that both S,rnM f^nd ^/fl.3£ individually approach 
zero as r 11/2, although their ratio in Eq. (|4.1ip 
approaches a constant. 

Illustrative values in the LNN limit for r = 4 are 



An 



,2e\r=i 



and 



m,3£k=4 



27r(-l + V5 
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1.396. 



0.8629. 



(4.12) 



so that for this value, r = 4, in addition to the ratio 
^.mad^iB^M = 2/3, we have 



r=4 



-1 + V5 



and 



£,m,3l 



= 0.68915 



0.6180. 



(4.14) 



(4.15) 



B. Value of at Minimum 

We calculate the following minimum values of P^^^ni as 
a function oi r £ Ir'. 



257r(ll - 2rf 
'34(13r-34)2 



(4.16) 



and 



37r(F3<> -I- GuVEm ) 



(4.17) 



where D^i and E^i were defined above in Eqs. (I2.26P 
and (|4.8p . The functions i^3^ and G31 are given in the 
appendix. 

For the illustrative value r — 4, 



ik.2i)r 



£1 
35 



= -0.1034 



and 



ik..3i)r 



27r(13- 5^5 ) 
35 



-0.04705 



(4.18) 



(4.19) 



C. dPi-^ni/di at ^/H,, 



At the two-loop level, we calculate 
2(11 -2r)2 



iJR,2e 



3(13r - 34) 



(4.20) 



This is clearly positive for r € Ir, approaching zero as r 
approaches the upper end of this interval at r = r'&iz- At 
the three-loop level, we find 



''^JR,3l 



4[-3(13r - 34)^3^ + i^3^VC^ ] 



^31 



(4.21) 



where Cs^ and D^^ were defined above in Eqs. (|3.9p and 
(I2.26p . and we define 



3t 



^21023 -f- 16557r - 3129r2 -|- 224r3 



(4.22) 



The first term in the numerator, — 3(13r — 34)C3£, is 
negative for r € Ir, but is smaller in magnitude than 
the second term, Ks^^/Csg. This shows analytically that 
(4.13) l3iiR.,3e > for r e /r. The positivity of P'^jR^ni for 
r g is obvious from the graph of (i^^ni- Since this 
function is continuous, is negative for < ^ < £,iR^ni, and 
has, generically, a simple zero at „^ > 0, it follows that 

l^'iRM > for r S Jr.. We have also calculated fi'^jnM 
alytically, but the expression is somewhat cumbersome, 
since it involves cube roots, so we do not list it. Illus- 
trative results of these ratios for the LNN limit and the 
typical value, r = 4, are 



and 



k, 



i,IR,3t\r=i 



X,IR,2l\r= 



4(-44-h 13V22 



0.2794 



(4.23) 



(4.24) 



This illustrates how the slope of the /3f function at the 
n-loop IR zero decreases as it is calculated to three-loop 
order, as compared with the 2-loop calculation. This is 
also shown by the explicit numerical results in the tables 
going up to four-loop order. 
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V. ANOMALOUS DIMENSION 7™ 

The anomalous dimension 7^ describes the scaling ol 
a fermion bihnear and the running of a dynamicahy gen- 
erated fermion mass in the phase with spontaneous chiral 
symmetry breaking. 7™ is defined as = din Zm/dt, 
where is the corresponding renormalization constant 
for the fermion bihnear operator. In the non-Abehan 
Coulomb phase (which is a conformal phase) , the IR zero 
of /3 is exact, although a calculation of it to a finite- 
order in perturbation theory is only approximate, and 
7m evaluated at this IR fixed point is exact. In the phase 
with spontaneous chiral symmetry breaking, where an IR 
fixed point, if it exists, is only approximate, 7™ is an ef- 
fective quantity describing the running of a dynamically 
generated fermion mass for the evolution of the theory 
near this approximate IRFP. As in [l^] , for notational 
simplicity we will often suppress the subscript m on 7^ 
where the meaning is clear. 

This anomalous dimension can be expressed as a series 
in a or equivalently, a: 



7 = 7" 



E 



cia 



E 



C£a 



(5.1) 



where £ denotes the loop order and q — Q/(47r)^ 
is the ^-loop series coefficient. The coefficient ci is 
scheme-independent while the ci with i > 2 are scheme- 
dependent. The ci have been calculated up to ^ = 4 in 
the AIS scheme The n-loop expression for 7 is de- 
noted jni and is given by Eq. (15. ip with 00 replaced by 
n as the upper limit on the summation over £. In Q we 
evaluated 7 to three- and four-loop order at the IR zero 
of /? calculated to the same order and showed that these 
higher-loop results were somewhat smaller than the two- 
loop evaluation. In [2^1 we discussed scheme transfor- 
mations at an IR fixed point and implications, including 
those for the scheme-dependence of 7. 

In the IR-conformal phase, unitarity implies a lower 
limit on the dimension of a spinless operator O, naniely, 
Do > (rf— 2)/2, where d is the spacetime dimension j27| . 
With our sign convention. 



In the IR phase with confinement and spontaneous chi- 
ral symmetry breaking, the dynamical mass generated 
for the fermions behaves, for Euclidean momentum large 
compared with the chiral-symmetry-breaking scale. A, 
as S(fc) cx A(A/fc)^^^ (up to a logarithmic factor) for 
fc >> A. Hence, the upper bound ()5.3p also applies in 
this phase. 

For our present analysis of the LNN limit, we reexpress 
7 in terms of x or equivalently, ^: 



(5.4) 



where 



and 



e=i 



ci = lim —T 

LNN M 



Ci = hm —J 

LNN Ni: 



(5.5) 



(5.6) 



Thus, similarly to the relation between a and q. 



Ce 



ci 



(47r) 



(5.7) 



The n-loop expression for 7 is given by the right-hand 
side of Eq. (j5.4l) with the sum running from £ = 1 to 
£ = n rather than £ = 00. 
In the LNN limit we find 



Cl 



C2 



203 
T2" 



(5.8) 



(5.9) 



D^,p — 3 



7 



(5.2) 



so that in this IR-conformal phase, 7 is bounded above 
as 



C3 



11413 



1177 



12C(3) 



35 
— 7 
27 



108 V 54 
105.676- 36. 221r- 1.296r2 (5.10) 



7 < 2 



(5.3) 



and 



460151 23816 899 , 83 , /1157 889 , 16 A,,,, 

C4 = r + + r + 2Gr^ H C(3) 

576 81 162 81 V 9 3 9 J^^ ' 

+ r(^66-12r)c(4)+ (- 220+ 160r)c(5) 

= 725.280 - 412.892r + 16.603r2 -|- 1.1123r3 , (5.11) 
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where the floating-point numerical results are given to 
the indicated accuracy. We list numerical values of the 



corresponding coefficients for 1 < ^ < 4 in Table IIVI 
In [9\ , the value of jni evaluated at a = am.ne was de- 



noted as , 
we define 



and, analogously, here, in the LNN limit. 



llR,. 



i—ilR.nt 



(5.12) 



At the two-loop level, in terms of the coefficients of f3 
and 7m, taking into account that 62 < for r d Ir, 



^ IR,2 



bi{ci\b2\ + C261) 
"2 



(5.13) 



Note that the sum of the i'-values of the products of co- 
efficients in the denominator of the expression for 7^^ 
is equal to the sum of the ^-values of the coefficients in 
the numerator. In this case, this sum is 4. Because of 
this homogeneity property and the relations in Eqs. (|2.6p 
and (|5.7p , it follows that 7^^ 3, has the same form as Eq. 
(|5.13p if one replaces each hi by hi and each ci by q. 
Exphcitly [9], 



(11 - 2r)(1009 - 158r + AQr'^ 



12(13r- 34)2 



(5.14) 



This 7^^ 2^ decreases monotonically as r increases from 
= 'rb2z = 34/13 at the lower end of the interval Ir to 
zero at r = rbi^ = 11/2 at the upper end of this interval. 
Because of the upper bound (|5.3p . if the value of 7 calcu- 
lated via this truncated pcrturbative expansion is greater 
than 2, then it is unphysical. This applies, in particular, 
to the divergence in 7^^ 3^ at r = rij2z- We find that in 
the relevant interval /r, ^,^21 exceeds 2 as r decreases 
below the value 3.569. Thus, we cannot use the formula 
(|5.14|) for r < 3.569, and it is subject to large corrections 
unless r is substantially above this value. For reference, 
we find that 7/i?,2£ exceeds 1 as r decreases below the 
value r = 3.879. 

At the three-loop level, 



(5.15) 



liR.ae = 2;(ci + C2X + CaX^) 



Substituting xm^^i — S,^^ 3^ / (47r) from Eq. (|3.6p and 
using the property that 62 < for r e Ir and, as dis- 
cussed above, the property that 63 < in a scheme that 
preserves the existence of the IR zero in the (scheme- 
independent) f3^^2£ at the three-loop level, we can write 
liR,M in terms of positive quantities as 



"flR,3 



461I63 



2C1I63P + C2I&2II63I 



• C3{bl + 261I63I) + (-C2I63I + C3\h2\)Jbl + 46i|5; 



(5.16) 



Thus, jjj^ 3j has the same type of homogeneity property 
as a function of the hi and C£ coefficients as we discussed 
for 7jjj 2f ■ Hence, 7jh,3<! is identically expressed by replac- 
ing each hi by hi and each q by C£ in Eq. ()5.16|) . Inserting 
the explicit expressions for the hi and C£ then yields the 
explicit result as a function of r. For the present work 
we have also calculated the four-loop anomalous dimen- 
sions evaluated at the IR zero of /S^ evaluated to the same 
order, 7,^ ^,. 

We list numerical values of these quantities in Table [V] 
We find that, as was the case with 7^-^ 1 liR.ae decreases 
monotonically as r increases from = rb2z at the lower 
end of the interval Ir to zero at r = r^i^ at the upper 
end of this interval. At r = ri,2z, Imai = 2.680, and 
hence the value of 3^ is evidently too large for the 
(perturbative) calculation of this anomalous dimension 
to be reliable. Jij^^i decreases through the value 2 as 
r increases through the value 2.73, and Jj^ ^i decreases 
further through the value 1 as r increases through 3.09. 
Just as was true of particular values of Nc and Nj Q, 



here, for a given value of r, 

"fiR.M < llR.2i > (5-17) 

with equality only at r = niz, where 7^^ 3^ = liR,2e = 0- 

VI. APPROACH TO THE LNN LIMIT 

One of the motivations for the present work on calcu- 
lations of the IR zero, the various structural properties 
of /?, and the anomalous dimension 7 evaluated at the 
IR zero in the LNN limit is that the results provide an 
understanding of (i) the approximate universality that is 
exhibited by calculations of these quantities for theories 
with different values of Nc and Nf such that the respec- 
tive values of r are the same or nearly the same, and (ii) 
the fact that this approximate universality occurs even 
for moderate values of Nc and Nf. As we have shown 
above, the quantities 

• aiR^niNc 



11 



(A 



da 



• 7n 



are finite in the LNN limit, and we will next show 
that they approach their respective LNN limiting values, 
namely 

• £,iR,ne 

min 



rather rapidly as Nc increases (with r — Nf/Nc fixed). 
The reason for this is that the largest correction terms to 
the respective LNN limits for these quantities are of or- 
der and hence are rather strongly suppressed even 
for moderate values of Nc and Nf . In turn, this is a con- 
sequence of the iVc-dependence of the relevant group in- 
variants for fermions in the fundamental representation. 



in particular, the Casimir invariant C/ = {N^ — l)/(2iVc) 



m p ; 



We first exhibit the structure of the correction terms 
for the coefficients in the beta function and anomalous di- 
mension 7. Multiplying by the appropriate inverse pow- 
ers of Nc to obtain finite results in the LNN limit, we thus 
consider bg/Nc and ci/Nc- These are polynomials in r 
and in the variable l/Nc- We find the following explicit 
results for correction terms: 



bi_ 

Nc 



= 61 



r 



(6.1) 



(6.2) 



llr(17- 2r) 
367V2 



4m 



(6.3) 



and 



1 



40 



58583r 2477H 



+ 



2341 
216 



1944 
623r 4 



243 



243 



352 



548r 64r" 



C(3) 



54 



(ll + 61r)C(3) 



23r 

STY? 



(6.4) 



For the ratios ce/Nc we find 



and 



C2 



Nc 



C2 ■ 



Cl 



3 



-53 + 5r 
3N^ 



1 /129 23r 

mKl, 2~ 



216 



12C(3)r - 



Am 



26309 899 



27 



129 

sm 



(6.5) 



(6.6) 



35 
27' 



(6.7) 



The corresponding result for C4 is given in the appendix. 
As noted, these expansions exhibit the feature that the 
largest subleading term is smaller than the leading term 
by a factor of 1/N^ in the LNN limit. 



Similarly, the approach of the product aiii^2iNc to its 
LNN limit, £,ir^21, has the form 



aiR,2eNc 



47r(ll - 2r) 127rr(ll - 2r) 
13r - 34 (34 - 13r)^N^ 



+ 



im) ' 



(6.8) 

where the first term on the right-hand side is the LNN 
limit, S,iR^2e, given in Eq. p.2p . Similarly, the two- loop 
anomalous dimension jm , evaluated at the IR zero of the 
two-loop /3 function, has the expansion 



(11 - 2r)(1009 - 158r + 40r2) 



12(13r- 34)2 
(11 - 2r)(18836 - 5331r + 648r2 - UOr^) 
(13r - 34)3iV2 



(6.9) 



where the first term on the right-hand side is the LNN 
limit, given in Eq. (j5.14p . Similar results hold for the 
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expressions calculated to higher-loop orders. Analogous 
results concerning correction terms to LNN limits apply 
to the theory with J\f — 1 supersymmetry, as will be clear 
from our discussion below. 

We next give some numerical results illustrating the 
rapidity of approach to the LNN limit. For definiteness, 
we consider two cases: Nc ~ 3, Nf = 12, so r = 4, and 
Nf, = 3, Nf = 15, so r = 5. For the first case, N^. ~ 3, 
Nf = 12, the two-loop, three-loop, and four-loop values 
of the IR zero of (3 are 

aiR,2e = 0.7540, ajR.si = 0.4349, = 0.4704 . 

(6.10) 

(In Table III of [9| these were listed these to three sig- 
nificant figures; here we list them to higher accuracy to 
compare with the estimates from the LNN limit.) To 
compute the LNN approximations to aiR,^nf. calculated 
for a specific Nc and Nf = riVc, we take the correspond- 
ing values for S^iR.ni from Table |llT] and divide by TV^, 
obtaining the result 



IR,ni 



for fixed r 



(6.11) 



where the subscript L indicates the LNN origin of the 
estimate. These approximations become progressively 
more accurate as Nc gets large for fixed r, but as we will 
show, they are already quite close to the actual values in 
Eq. dmni). We obtain 

aiR,2e.L = 0.6982, aiR,3i,L = 0.4713, a/K,«,L = 0.4497 

(6.12) 

Hence, for this Nc = 3, Nf — 12 case, we find 



Performing the corresponding analysis for Nc = 3, Nf 
15, we obtain 



lIR,2t 
1IR,21X 



= 0.943, 



0.929, -^^^^^= 0.929. 

1IR,U,L 



3.17) 

These numerical comparisons show that even for a mod- 
erate value of Nc such as iV^ = 3, the values of these IR 
zeros of /3q,„,£ and anomalous dimensions jiR^ne are close 
to the approximations that one obtains from the LNN 
limit. The agreement becomes better as Nc increases. 

From our analysis, it follows that the approach to the 
LNN limit is of the same form for other structural prop- 
erties describing the UV to IR evolution of the theory; 
that is, the first subleading correction term to the LNN 
limit is suppressed by l/Nc- We show this explicitly for 
additional two-loop quantities. The result for am,2i fol- 
lows immediately from Eqs. (|4.5p . (I4.6p . and (|6.8p . For 
{PiR,2i)mim the approach to the LNN limit has the form 



{PlR,2i)minNc — {/3^,2i)r. 



Mm 



33(13r - 34)3^2 ■ " Va4; ' 

(6.18) 

where {(3^^2e)min was given in Eq. (j4.16p . For the deriva- 
tive of (3 at the n-loop IR zero of /?, the approach to the 
LNN limit has the form 



2£ 



21 



da 

where dl3^^2e/d£, 

Km . 



i—^IR,'. 



2r(ll-2r)2 / 1 \ 
(13r - 34)2iV2 ^ \N^J 
(6.19) 



— I3'^ 21 given in Eq. 



aiR^2i 



1.080, 



aiR.Si 



1.042, 



1.046, 



Q:/_R,2£,L ClIR.Zl,L OliR^M.L 

(6.13) 

Carrying out the corresponding analysis for Nc — 3, 
Nf ~ 15, we calculate 



Q-IR,2t 



= 1.057, 



.^£^^ 1.046, ^1^ = 1.046 



OlIR,2l,L a:iR,u,L aiR^M.L 

(6.14) 

We next perform the corresponding numerical compar- 
isons for "fiR^ni at the two-loop, three-loop, and four-loop 
levels. For iV^ = 3 and Nf = 12, 

liR.2e = 0.7728, 7/^,3^ = 0.31175 7/^,4, = 0.2533 . 

(6.15) 

(In Table VI of 19| these were listed to three significant 
figures; here we list them to higher accuracy to compare 
with the estimates from the LNN limit.) The LNN values 
of 

^iR,ni are displayed in Table IVl Comparing these with 
the values in Eq. (pTTSt . we find, for this Nc = 3,Nf = 12 
case, the ratios 



VII. CONTRAST WITH FERMIONS IN 
HIGHER REPRESENTATIONS 

Gauge theories with fermions in higher-dimensional 
representations, in part icular, two-index representations, 
are also of interest [28| . We comment bricfiy here on the 
adjoint, and symmetric and antisymmetric rank-2 tensor 
representations with Young tableaux I I I and Rather 
than an LNN limit, in these cases, one takes Nf equal 
to a (non-negative, integer) constant as Nc — 00 to get 
a finite limit. For the adjoint representation, bi/Nc is a 
constant for ^ — 1, 2, 3, while b/i/Nc is equal to a con- 
stant plus a 1/Nc correction term (where hg with £ — 3, 4 
are calculated in the MS scheme). For the S2 ( m ) and 
A2 (Q) representations, symbolized together as T2, the 
corrections to the Nc 00 limit go like 1/A^c instead of 
l/Nc- For example. 



Nc 



11 - 2A^ 



2Nf_ 

Nc 



for T2 , 



(7.1) 



'^^^■^^ - 0.985, ^^^'^^ 



lIR,2e,L 



1IRM,L 



0.913, ^^^.4^ ^ 



1IR,4e,L 



3.16) 



where the upper (lower) sign applies to S2 (A2). 

These differences are reflected in the large- A''^, correc- 
tions to the n-loop expressions for the IR zero of /3. We 
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illustrate this at the two-loop level. For the adjoint rep- 
resentation, Nf^biz = 11/4 and Nf^i,2z — 17/16, so that 
there is only a single integer value of Nf for which the 
theory is asymptotically free and has an IR zero in P2i, 
namely Nf = 2. One has 



27r(ll-4iVf) 



(7.2) 



(independent of Nc) so the right-hand side is equal to 
27r/5 for the value Nf = 2, independent of A^^- 

For the S2 and A2 representations (denoted T2 again), 
for large N^ 



2n{ll-2Nf) 6TTNf{-47 + 2Nf) 



87V/ - 17 



(87V/ ~ llfN, 



o(-K\ for R = T2 



(7.3) 



so that in these S2 and A2 cases, the leading correction 
to the Nf, ^ oo result goes like l/N^- 



VIII. SUPERSYMMETRIC GAUGE THEORY 
A. /Jj s Function and IR Zeros 

Here we consider the LNN limit of an asymptotically 
free, vectorial gauge theory with J\f = 1 supersymmetry, 
gauge group G = S\J{Nc), and a chiral superfield con- 
tent <I>i, $i, i — l,...,Af/, in the □, □ representations, 
respectively. One of the appeals of this theory is that 
a number of exact results on the infrared properties of 
the theory are known [2^, [sO] , so one can compare per- 
turbative predictions with these exact results. This was 
done for general G and various representations^ R for 
the Nf pairs of chiral superfields $i, in [ill, Our 
discussion of the LNN limit here extends the previous 
results in [Til, [13 (see also ^,^). 

The /3 function of this theory will be denoted /J^^s and 
has the expansion (|2.2p with bg and bg replaced by be,s 
and hi s- Here and below we use the subscript s, standing 
for "supersymmetric" , to avoid confusion with the cor- 
responding quantities discussed above in the nonsuper- 
symmetric theory. Thus, 



be.s = lim , 6f.., = lim 



LNN 7V5 



LNN Ni 



(8.1) 



We denote the n-loop /3^_s function as P^^sM- The 
scheme-independent coefficients bi^s and 62. s, were com- 
puted in ,325 and [11], and 63^^ was computed in [sj] in 
the dimensional reduction {DR) scheme [35| . 
In the LNN limit, 



bi,s = 3 



2(3 - 2r) , 



(8.2) 



(8.3) 



and, in the DR scheme, 

= 21 - 21r 



4r^ 



(8.4) 



From these values of the coefficients of /3j_s, it follows 
that 



and 



n2z 



(8.5) 



(8.6) 



Asymptotic freedom requires r < r^i^, and, in this range, 
the two- loop 13 function has an IR zero for r > ri,2z, so 



the interval Ir « is 



< r < 3 



(8.7) 



This IR zero occurs at the value 



21,! 



hi^s _ 47r6i^, 

&2.S h2.s 

27r(3- r) 



2r - 3 

The coefficient hz,s vanishes at two values, 
_ 21 ± \/T05 

^fc3^,s,(l,2) — 7, ■ 



(8.9) 



Numerically, rf,3z,s,i — 1.344 and rt,3z,s.2 — 3.906, so that 
rb?.z,s,l < n2z , ?'632,s,2 > niz ■ (8.10) 

Since 63 ^ < for r^s^ s 1 < r < ri,-i,z.s,2, it follows that, 
in the DR scheme. 



< V r eir 



(8.11) 



We list numerical values of the bi,s in Table IVIl 

The three-loop /? function formally vanishes at two 
points away from the origin, at the zeros of bi g + 62, + 
63, sC^, namely 



1 



2b 



3,s 



± 



bis 



'3,s 



5.12) 



By the same type of argument that was given above for 
the nonsupersymmetric theory, one may argue that the 
inequality (18. lip applies more generally than just in the 
DR scheme. In the present context, this argument is that 
if the theory has an IR zero in the (scheme-independent) 
two-loop /? function, it is reasonable to require that a 
scheme should preserve the existence of this IR zero at 
higher-loop level, this requires that 63 < Vr € I^-^s- 
This follows since 62,3 at the lower end of the inter- 
val Ir.s^ so unless &3 s < in this interval, the quantity 
^2 s — 46i_s&3,s would become negative and the value of 
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the zeros of /3 at the three-loop level in Eq. (j8.12p would 
be complex. Given that &3 < for r € Ir,s, one may 
rewrite Eq. (|8.12p in terms of positive quantities and 
pick out the relevant IR zero of (3^.s,3t as 



IB.,3l,s 



Explicitly, 



216: 



3,s| 



02, s| 



52 +46i,,|53,s| 



4tt 



^IR,U,. 



(2r - 3) + 



where 



and 



= -54 + 72r - 2<dr'^ + Ar^ 



-21 + 21r-4r^ 



.13) 



(8.14) 



.15) 



.16) 



Note that Ds > for ri,3z,s,i < r < rb3z.,s,2 and hence 
for all r g Ir^s- Furthermore, Cs is positive-definite for 
r e Ir^s (the zeros of Cs occur at r = 1.338 and r — 
2.956 ± 1.163«). By the same reasoning as was given 
before in and above, we have the inequality 



i.n) 



We list numerical values oi£,iB.,ni-s i^i Table IVlIl and show 
a plot of P(^ne,s for n — 2 and n = 3 loops, evaluated at 
an illustrative value r in Ir^s, namely, r = 2.5 in Fig. [51 



B. Other Structural Features of P^^nt.s 

We also briefly discuss some structural features of f3^,s- 
The two-loop /3 function, (3^^s,2i, reaches a minimum on 
the interval Ir s at the value 



.2f,, 



47r(3- 



2r- 3 



(8.18) 



The ratio of this position of the minimum in the beta 
function relative to the position of the IR zero is the 
same as in the nonsupersymmetric theory, namely 



^IR,2l,s 3 

The value of P^,2i,s at the minimum is 



21,. 



87r(3- 



?.19) 



(8.20) 



27(2r-3)2 ■ 

The derivative df3^j[i^2i.s evaluated at ^ = £.iR,2i.s, is 
_ (3 - r)2 



i,IR,2i,s 



2r-3 



(8.21) 



Corresponding expressions can be given at the three-loop 
level, but we proceed now to analyze a quantity of con- 
siderable interest, namely the anomalous dimension of 



C. Anomalous Dimension 

We next consider the LNN limit of the anomalous di- 
mension of the (gauge-invariant) quadratic chiral super- 
field product, $$, denoted 7s = jm,s- This is given by 
Eq. (15.11) with the replacements ci — C£^s and cg — ^ C£ 
and similarly for the n-loop expression, jrn,,s,ni- From 
the known results for ci.s and, in the DR scheme, C2,s 
and C3 s, we find, in the LNN limit. 



ci s = 2 



C2,, 



2(2 - r) 



and 



C3,; 



10-6r[l + 2C(3)] -2r^ 



..22) 



5.23) 



^.24) 



Values of the corresponding ci^s are listed as a function 
of r in TableEnH 

The two-loop anomalous dimension, evaluated at the 
two-loop IR zero of is [ll| 

r(r — 1)(3 — r) 

The quantity 7^^ ^ is a monotonically decreasing func- 
tion of r in the interval Ir,s- It exceeds the upper bound 
of 1 from conformal symmetry if r < 2, and hence, as 
was noted in [ll|, in the interval 3/2 < r < 2, the per- 
turbative two-loop calculation that yields ^jf, 2e,B gives an 
unphysical result and is unreliable. One can apply this 
upper bound because one knows from exact results [s^l 
that for 3/2 < r < 3, the theory flows in the infrared to 
a conformal, non-Abelian Coulomb phase. 

The three-loop anomalous dimension, evaluated at the 
three- loop IR zero of is 



li 



(8.26) 



where Cs and Ds were defined above in Eqs. (I8.15P and 
(I8.16|) . and As and Bs are given in the appendix. In 
contrast to 7^^ 3 : Iir u s is not a monotonic function 
of r. It reaches a maximum of approximately 0.1376 at 
r — 2.474 and vanishes not just at r = rbiz = 3, but also 
at r = 2.1794 (to the indicated number of significant fig- 
ures). These features were evident for the specific values 
of Nc considered in [ll[; here we have shown how this 
occurs after the LNN limit is taken. In Table |TX] we list 
values of 7^^ , for a range of r values in Ir.s- We note 
that for 1.5 < r < 2.0, Jj^f ^i 3 exceeds the upper bound 
of 1 and hence is unphysical; we do not include entries 
for these values of r. 



IX. DISCUSSION AND CONCLUSIONS 

In summary, we have studied higher-loop corrections 
to the UV to IR evolution of an asymptotically free vecto- 
rial SU(7Vc) gauge theory with Nf fermions in the funda- 
mental representation, in the 't Hooft-Veneziano (LNN) 
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limit Nc oo and Nf — oo with r = Nf/Nc fixed 
and ^(^) — a{^)Nc, a function independent of in this 
hmit. We have defined a beta function, /3^, that is finite 
in this LNN Hmit and have analyzed its properties for 
the interval of r in which /3{,2£ has an IR zero. We have 
given analytic and numerical results for the LNN lim- 
iting quantities ^/_r,„£, ^m.ni^ il3^.ni)min, P'^.iR.ni^ and 
liR.ni as functions of r. We have argued that a rea- 
sonable scheme should preserve at higher loops the IR 
zero that is present at the (scheme-independent) two- 
loop order in /3, and that this implies that &3 < for 
r £ Ir- In turn, this implies that S,ir^u < iiR,2i- Cal- 
culating with the MS scheme, we find that in the part 
of the interval Ir where the perturbative calculations are 
reliable, the change in the IR zero is smaller in magni- 
tude going from three-loop to four-loop order, as com- 
pared with the shift from two-loop to three-loop order. 
This is in agreement with one's expectation, that in- 
sofar as perturbative methods are trustworthy, when a 
quantity is calculated to higher orders, the successive 
changes should become smaller. Further, we find that 
for the range of r S /r where the three-loop anomalous 
dimension is reliably calculable, ^ir^u < liR,2i- These 
higher-loop calculations allow one to extend the analysis 
of the IR zero of the /3 function, and corresponding eval- 
uations of 7m to smaller values of r and thus stronger 
couplings than is possible with the two-loop result. We 
have analyzed the correction terms to the LNN-limit val- 
ues of a number of quantities, including hi/N^, cg/Nf^, 

OtlR^neNc, am,nl^c, {l3a,ni)ininNc, dl3a,ni/d'^\aiR.nn and 



liR.nl, and have shown that these correction terms are 
suppressed by 1/iV^. This provides an understanding of 
the approximate universality that is exhibited in calcula- 
tions of these quantities for different values of Nc and Nf 
with similar or identical values of r, even for moderate 
values of Nc and N f. A corresponding analysis was also 
given of a vectorial gauge theory with M = 1 supersym- 
metry and Nf chiral superfields transforming according 
to the fundamental and conjugate fundamental represen- 
tation of SU(A^c)- Thus, in addition to being of interest 
in its own right, the LNN limit is useful in understand- 
ing common features of the UV to IR evolution of various 
theories with different values of Nc and N f. 
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X. APPENDIX 

In this appendix we list some of the more lengthy ex- 
pressions that are used in the text. The functions F^f^ 
and Gzt that enter in Eq. (|4.17|) for {fi^,zi)min are 



Fu = -10985980784+ 17408705952r - 10177907376r-' + 2883132208r-' - 468256107r^ + 42131712r^ - 1605632r'' 



and 



Gu = 41737992 - 48660252r -I- 18696078^2 - 2733297r3 -|- 139776r'' . 



(10.1) 
(10.2) 



For the coefficient C4 in Eq. (|5.1|) (with fermions in the fundamental representation, and in the MS scheme), using 
[ii, we find 



C4 = C4 





1 




21947 


126689 


C4 + 


iVf 




24 + 


324 ' 


1 


'108359 


9143 


38 2 
H 


Nj 


288 


r 

108 


27 


1 


5783 


37 


(89 + 1 






24 





162 



83 , / 221 379 16 ^\^,^, ,,,, 
+ + ( - — + — - )C(3) - 60rC(5) 



1 



1261 
64 



- 42C(3) 



(10.3) 



The functions Ag and Bg that enter in the expression for 7^-^ 3^ ^ are 
Ag ^ 918 - 2916r + AlAQr^ - 3322r^ + 1532r* - 378r^ + 40r*^ + r ^2754 - 5508r + 3942r2 - 1212r^ + 144r^ ) C(3) (10.4) 



and 



Bg^M- 75r + 17r2 + + br'^ - Ar^ + r[lQ2 - 216r + lQ2r'^ - 2Ar^^ C(3) • 



(10.5) 
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TABLE I: Values of the coefficients for 1 < ^ < 4 as functions 
of r for < r < ri,!^. Notation ae-n means a x 10~" here and in 
the other tables. 



r 


bi 


62 


63 


64 


0.0 


0.2918 


0.7177e-l 


0.2666e-l 


0.1265e-l 


0.5 


0.2653 


0.5805e-l 


0.1895e-l 


0.8063e-2 


1.0 


0.2387 


0.4433e-l 


0.1176e-l 


0.4429e-2 


1.5 


0.2122 


0.3061e-l 


0.5091e-2 


0.1766e-2 


2.0 


0.1857 


0.1689e-l 


-0.1055e-2 


0.9062e-4 


2.5 


0.1592 


0.3166e-2 


-0.6677e-2 


-0.5814C-3 


3.0 


0.1326 


-0.1055e-l 


-0.1178e-l 


-0.2340e-3 


3.5 


0.1061 


-0.2427e-l 


-0.1635e-l 


0.1149e-2 


4.0 


0. 7958c- 1 


-0.3800C-1 


-0.2041C-1 


0.3584e-2 


4.5 


0.5305e-l 


-0.5172C-1 


-0.2394C-1 


0.7086C-2 


5.0 


0.2653e-l 


-0.6544e-l 


-0.2695e-l 


0.1167e-l 


5.5 





-0.7916e-l 


-0.2944e-l 


0.1736e-l 



TABLE II: Values of the with 1 < ^ < 4 at special values of r, 
including r = and at the lower and upper ends of the interval Ir , 
r = n2z = 34/13, and r = = 11/2. 



£ 


ibi)r=0 


ibe)r = rt2;. 




1 


0.2918 


0.1530 





2 


0.7177e-l 





-0.7916e-l 


3 


0.26660-1 


-0.7900C-2 


-0.2944e-2 


4 


0.1265e-l 


-0.5923e-3 


0.1736e-l 
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TABLE III: Values of the IR zero (iR^ni of /9£,n< function for 
n = 2, 3, 4 and r e Ir- 



r 


£,IR,2t 




£,IRAI 


2.8 


28.274 


3.573 


3.323 


3.0 


12.566 


2.938 


2.868 


3.2 


7.606 


2.458 


2.494 


3.4 


5.174 


2.076 


2.168 


3.6 


3.731 


1.759 


1.873 


3.8 


2.774 


1.489 


1.601 


4.0 


2.095 


1.252 


1.349 


4.2 


1.586 


1.041 


1.115 


4.4 


1.192 


0.8490 


0.9003 


4.6 


0.8767 


0.6725 


0.7038 


4.8 


0.6195 


0.5083 


0.5244 


5.0 


0.4054 


0.3538 


0.3603 


5.2 


0.2244 


0.2074 


0.2089 


5.4 


0.06943 


0.06769 


0.06775 



TABLE IV: Values of the q coefficients with 1 < ^ < 4 as func- 
tions of r for < r < rji^. Here ci = 3/(47r) = 0.2387, indepen- 
dent of r. 



r 


C2 


C3 


C4 








0.1071 




0.5325e-l 


0.2908e-l 





5 


0.1018 




0.4396e-l 


0.2098e-l 


1 





0.9657e- 


1 


0.3434e-l 


0.1324e-l 


1 


5 


0.9129e- 


1 


0. 2440c- 1 


0.5897e-2 


2 





0.8602O- 


1 


0.1413C-1 


-O.lOlOe-2 


2 


5 


0.8074e- 


1 


0.3538C-2 


-0.7451e-2 


3 





0. 75460- 


1 


-0.7384C-2 


-0.1339e-l 


3 


5 


0.7019O- 


1 


-0.1863C-1 


-0.1880e-l 


4 





0.6491c- 


1 


-0.3021C-1 


-0.2364e-l 


4 


5 


0.5963c- 


1 


-0.4211C-1 


-0.2788C-1 


5 





0.5435e- 


1 


-0.5434e-l 


-0.3148e-l 


5 


5 


0.4908e- 


1 


-0.6690e-l 


-0.3442e-l 



FIG. 1: Plot of the n-loop beta function, (^^^n£, as a function of 
5, for the illustrative value r = 4. Prom bottom to top, the curves 
represent /3j,2^i ^■iid /3j,3f> respectively. 
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TABLE V: Values of the n-loop anomalous dimension, 7„f , eval- 
uated at the n-loop IR zero of /3j and denoted , as in Eq. 
I I5.I2I I. where n = 2, 3, 4, for r £ I^. In the entries marked u, the 
n-loop perturbative value of ■yjji „i is larger than the upper bound 
of 2 and hence is unphysical (u). 



r 


7lH 21 






2.8 


u 


1.708 


0.1902 


3.0 


u 


1.165 


0.2254 


3.2 


u 


0.8540 


0.2637 


3.4 


u 


0.6563 


0.2933 


3.6 


1.853 


0.5201 


0.3083 


3.8 


1.178 


0.4197 


0.3061 


4.0 


0.7847 


0.3414 


0.2877 


4.2 


0.5366 


0.2771 


0.2664 


4.4 


0.3707 


0.2221 


0.2173 


4.6 


0.2543 


0.1735 


0.1745 


4.8 


0.1696 


0.1294 


0.1313 


5.0 


0.1057 


0.08886 


0.08999 


5.2 


0.05620 


0.05123 


0.05156 


5.4 


0.01682 


0.01637 


0.01638 



TABLE VI: Values of the bi g coefficients for 1 < £ < 3 as functions 
of r for < r < rbi2,s in the supersymmetric theory. 



r 


bi,s 


62,5 


&3,s 


0.0 


0.2387 


0.3800e-l 


0.1058e-l 


0.5 


0.1989 


0.2533e-l 


0.5795e-2 


1.0 


0.1592 


0.1267e-l 


0.20166e-2 


1.5 


0.1194 





-0.7559e-3 


2.0 


0.7958e-l 


-0.1267e-l 


-0.2520e-2 


2.5 


0.3979e-l 


-0.2533e-l 


-0.3276e-2 


3.0 





-0. 3800c- 1 


-0.3024e-2 



TABLE VII: Values of the IR, zero 5/fl,„£,s of /3^,„£,s function for 
n = 2 and n = 3 loops and r £ Ir,s- 



r 


£,ir;21,s 




1.8 


12.566 


5.331 


1.9 


8.639 


4.381 


2.0 


6.283 


3.643 


2.1 


4.712 


3.040 


2.2 


3.590 


2.529 


2.3 


2.749 


2.085 


2.4 


2.094 


1.692 


2.5 


1.571 


1.339 


2.6 


1.142 


1.019 


2.7 


0.7854 


0.7279 


2.8 


0.4833 


0.4623 


2.9 


0.2244 


0.2201 



-0.1 J 



FIG. 2: Plot of the n-loop beta function, P^^ril.s, as a function 
of for the illustrative value r = 2.5 in the supersymmetric the- 
ory. Prom bottom to top, the curves represent /3j,2£,s ^^'^ Pi,3l,s^ 
respectively. 



TABLE VIII: Values of the ^ coefficients for 1 = 1 and £ = 3, 
as functions of r for < r < rfji^^s in the supersymmetric theory. 
Here cx^^ = l/(27r) = 0.1592. independent of r. 









r 






0.0 


0.2533e-l 


0.5039e-2 


0.5 


0.1900e-l 


-0.3590e-3 


1.0 


0. 12670-1 


-0.6261e-2 


1.5 


0.6333c- 2 


-0.1267C-1 


2.0 





-0.1958e-l 


2.5 


-0.6333e-2 


-0.2699e-l 


3.0 


-0.1267e-l 


-0.3491e-l 
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TABLE IX: Values of the n-loop anomalous dimension, 7„^^3 , eval- 
uated at the n-loop IR zero of /9j_s and denoted 7jjj ^ , where 
n = 2 and n = 3, for r 6 Ir^s- For 1.5 < r < 2.0, 7jjj 2< a ^ ^ ''■"'^ 
hence is unphysical (u). 



T 




"fiR,:u,s 


1 o 
i.O 


u 


— i.Dl / 


1 n 

i.y 


u 


n OAK Q 


2.0 


1 


-0.3667 


2.1 


0.7219 


-1183 


2.2 


0.5388 


0.2267e-l 


2.3 


0.4088 


0.9809e-l 


2.4 


0.3111 


0.1314 


2.6 


0.2344 


0.1370 


2.7 


0.1719 


0.9955C-1 


2.8 


0.7456e-l 


0.6829e-l 


2.9 


0.3514e-l 


0.3412e-l 



